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We investigate the time evolution of the temperature and entropy of a gravitationally collapsing 
cylinder, represented by an infinitely thin domain wall, as seen by an asymptotic observer. Previous 
work has shown that the entropy of a spherically symmetric collapsing domain approaches a constant, 
and we follow this procedure using a (3+1) BTZ metric to see if a different topology will yield 
different results. We do this by coupling a scalar field to the background of the domain wall and 
analyzing the spectrum of radiation as a function of time. We find that the spectrum is quasi- 
thermal, with the degree of thermality increasing as the domain wall approaches the horizon. The 
thermal distribution allows for the determination of the temperature as a function of time, and 
we find that the late time temperature is very close to the Hawking temperature and that it also 
exhibits the proper scaling with the mass. From the temperature we find the entropy. Since the 
collapsing domain wall is what forms a black hole, we can compare the results to those of the 
standard entropy-area relation. We find that the entropy does in fact approach a constant that is 
close to the Hawking entropy. However, the time dependence of the entropy shows that the entropy 
decreases with time, indicating that a (3+1) BTZ domain wall will not collapse spontaneously. 



I. INTRODUCTION 

It is well known, based primarily on the work of Beken- 
stein, Gibbons, and Hawking, that the entropy of a black 
hole is proportional to its area and that even though 
supposedly nothing can escape from within a black hole, 
quantum fluctuations near the event horizon will produce 
a spectrum of radiation that is thermal [IH3] . Since then, 
much work has been done to reproduce this result with 
theories of quantum gravity, but most of these theories 
do not analyze the time evolution of the system. The 
conventional process to determine the entropy uses the 
Bogolyubov method to first determine the temperature 
and from there find the entropy. This, however, only 
utilizes the initial and final states of the system, and 
therefore there is no knowledge of the time dependence. 
Recently, Vachaspati and Stojkovic developed a quantum 
treatment to determine the quantum radiation given off 
during gravitational collapse in a time dependent man- 
ner, then Greenwood expanded on this to determine the 
time dependence of the entropy @H1] ■ These papers used 
a spherically symmetric collapsing domain wall for their 
analysis, and we wish to augment their work by seeing if a 
different topology will yield significantly different results. 
Specifically, we wish to determine the time dependence 
of the entropy of a gravitationally collapsing cylinder, 
which we will represent with an infinitely thin domain 
wall. Ref.fTT] used a (3+1) BTZ metric to determine the 
equation of motion of the collapsing cylindrical domain 
wall, and we will use that to complete the analysis of 
its thermodynamic properties. To do this, we will first 
determine the wavefunctional of a scalar field coupled 
to the background of the collapsing shell. Then, using 
the t=0 wavefunctional as a basis, we will determine the 
occupation number as a function of frequency. This in 
turn will allow for the determination of the thermody- 
namic quantity j3 and therefore the temperature. We 
then determine the late time temperature as a function 
of horizon radius. Using the thermodynamic definition 



of entropy, dS = dQ/T, we integrate to find entropy as 
a function of temperature, and since we know the tem- 
perature as a function of time we can find the entropy 
as a function of time. We compare our results against 
the standard Hawking results, and we comment on the 
results. 



II. RADIATION WAVEFUNCTIONAL 

First, we will consider the radiation given off by the 
cylindrical domain wall during gravitational collapse. 
This will be done by coupling a scalar field, which we 
will decompose as 



$ = ^a k {t)u k {r), 
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to the background of the collapsing shell. The exact form 
of Uk(r) will not be important to us. To find the modes 
Ofe(i), we will insert metrics (see Ref. |llj) 
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The total action will be written as the sum 
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is the metric outside the shell, and 
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In the region of interest, where R — > Rr and therefore 
/ — > 0, the kinetic term in Si n dominates over that in 
S out while the gradient term in Sj, n is subdominant to 
that in Sout- This yields the approximate action 
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The action ( 18 1 corresponds to the Hamiltonian 
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We can now find the wave function ip from 
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Since M and N are Hermitian matrices, it is possible to 
do a principle axis transformation to diagonalize them 
simultaneously. The exact form of them, however, is not 
important here. The Schrodingcr Equation for a single 
mode will then be 
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where m and K are the eigenvalues of M and N, respec- 
tively, and b is the eigenmode. 

This, however, only represents the Hamiltonian of the 
radiation, not the total system. The Hamiltonian of the 
system can be written as 



Htotal — H wa u + H ra d- 

Ref.[TT] showed that 
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where fj, = \J AR 2 H /3 — 2-koGRh and a is the tension in 
the wall. As R — >• / — > 0, and the first term on the 
right hand side will dominate (since n ~ J -3 ^ 2 ), yielding 
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where we have chosen the negative solution because the 
shell is collapsing. The Hamiltonian of the radiation is 
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The Hamiltonian of the entire system is therefore 
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Note that the momentum operators have different in- 
dices, where = —id/dR is the momentum operator 
conjugate to the shell's position and Hi, = —id/db is the 
momentum operator conjugate to the eigeinmode b. The 
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wavefunction of the entire system is therefore a function 
of 6, R 7 and t, which we will write as 
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(28) 



Using (|27|), the Functional Schrodinger equation can be 
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Solving this requires knowing the classical equation of 
motion for the distance of the domain wall, which we 
will take as approximately 
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where again the negative sign was chosen because the 
shell is collapsing. Since R(t) is only a function of time, 
we can rewrite ( 29 ) as 
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According to (30), however, f/R = —1, therefore 
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Rewriting this in the standard simple harmonic oscillator 
form yields 
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Here we have chosen r](t = 0) = 0. 

To proceed further, we need to incorporate the equa- 
tion of motion R(t) given by (see Ref.|llj) 



R(t) = R h + {Ro - R h )e 
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where we are assuming a negative cosmological constant. 
At early times, R(t) is approximately constant, so the 
initial vacuum state for the wavefunction tp is that of a 
simple harmonic oscillator given by 

tf(M = 0)=(^) 1/4 e-™- °" 2 / 2 . (37) 
At later times, 
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and the exact solution for the wavefunction is (see 
Ref.p]) 
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where p n denotes the derivative of p(rj) with respect to 
rj, and p(rf) is the solution to the equation 
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The phase a is given by 
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III. OCCUPATION NUMBER OF THE 
RADIATION 

Consider an observer with detectors that are designed 
to register particles for the scalar field (f> at early times. 
At late times, the observer will interpret each mode b 
in terms of the simple harmonic oscillator states, with 
final frequency u). The number of quanta in eigenmode 



b can be found by decomposing the wavefunction ( 39 ) 



into the simple harmonic oscillator states and evaluating 
the occupation number. The wavefunction ip written in 
terms of the simple harmonic basis, {f n }, at t = is 
given by 



where 
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which is the overlap of a Gaussian with the simple har- 
monic basis functions. The occupation number at eigen- 
frequency Co is given by the expectation value 



(46) 



After substitution, we find that the occupation number in 
the eigenmode b is given by (see Appendix B in Ref.|13|) 
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The plot in Fig.Q shows the occupation number as a 
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FIG. 1: Occupation number of the radiation as a function of 
frequency. From bottom to top, the three lines are for t=13, 
t=14, and t=15, respectively. 
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FIG. 3: /3 as a function of time for a collapsing 3+1 BTZ 
domain wall. 



function of frequency for three time slices. This is very 
similar to a typical Planck distribution. Therefore, we 
can find the temperature of the radiation by comparing 
the occupation number of each eigenmode b to that of 
the Plank distribution, 

N * = ( 48 ) 

If one plots ln(l + 1/N) as a function of u>, the slope 
will yield j3, as shown in Fig. It should be noted that 
this slope is conjugate to the coordinate r\. If we wish to 
obtain /3 as a function of observer time t, then we must 
rescale according to 
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FIG. 2: The straight lines are the best fit to the curve and 
the slopes represent f3. 

temperature as a function of time, as shown in Fig.Q. 
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FIG. 4: Temperature as a function of time for a collapsing 
3+1 BTZ domain wall. 



We will now compare the late time temperature to the 
Hawking temperature. According to |15j . 




The ratio of the late time temperature in the plot to 
the Hawking temperature for Rh = 1 is Tbtz/Th = 
1.33. The plot of the domain wall's late time temperature 
versus mass, as well as the best-fit curve, is shown in 
Fig. ([5j . The equation of the best-fit curve is 

Tbtz = 0.131724M 1/3 (51) 

We can therefore see that the late time temperature ex- 
hibits the proper scaling with the mass of the black hole, 
that is T ~ M 1 / 3 . 
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FIG. 5: Temperature as a function of mass for a BTZ domain 
wall. 



IV. ENTROPY 

The thermodynamic definition of entropy in terms of 
temperature is 



S = 



dQ 
T ' 



(52) 



Since changing the energy of the domain wall is the same 
as changing the mass, this can also be written as 



8 = 



dM 
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Therefore, if we can determine how the temperature of 
the domain wall depends on the mass, we will be able to 
find an expression for the entropy in terms of the tem- 
perature. And since we know the time evolution of the 
temperature, we can determine the time evolution of the 
entropy as well. 



We insert Eq.(51) into Eq.(53), and after integrating 



we find that the entropy as a function of temperature is 
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where 7 = 0.131724. 

The corresponding plot of entropy as a function of time 
is shown in Fig. ([6]). Notice that since the entropy is pro- 
portional to the square of the temperature, the entropy 
decreases along with the temperature. At late times, 
the domain wall reaches a static state and the entropy 
ceases to decrease. According to Ref [TB] . a BTZ black 
string does not evaporate once formed, and is thus sta- 
ble. However, our result of decreasing entropy seems to 
imply that a (3+1) BTZ domain wall would not collapse 
spontaneously under these conditions, so the black string 



would never be formed in this way. In fact, the assump- 
tion of a collapsing cylinder may be wrong entirely, mean- 
ing that perhaps it is the expanding solution which leads 
to an increase in entropy. 

Sbtz 




FIG. 6: Entropy as a function of time for a BTZ domain wall. 



V. CONCLUSION 

We investigated the time evolution of the temperature 
and entropy of a 3+1 BTZ domain wall (representing a 
cylinder) . This was done by coupling a scalar field to the 
background of the domain wall and evaluating the occu- 
pation number. From the occupation number we were 
able to determine /3 and therefore the temperature. We 
found that the temperature of the 3+1 BTZ domain wall 
decreased, then approached a constant, and that the late 
time temperature exhibited very good agreement with 
the Hawking temperature of a static black hole. It also 
produced scaling with mass that matched the literature. 
The entropy, however, exhibited very interesting behav- 
ior, as the entropy actually decreased over time. At this 
point we are unable to determine the exact cause of the 
decrease, however, because this particular metric con- 
tains both an anti-de Sitter cosmological constant and 
a different topology from previously studied cases. The 
interesting implication here, though, is that the collapse 
of a (3+1) BTZ domain wall does not happen sponta- 
neously. Future work could examine whether the expand- 
ing solution produces an increasing entropy. 
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